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ABSTRACT 

We have obtained initial spectroscopic observations and additional photometry of the newly discovered Pf, — 
94min 7-ray black-widow pulsar PSR J131 1—3430. The Keck spectra show a He-dominated, nearly H-free 
photosphere and a large radial-velocity amplitude of 609.5 ± 7.5 km s^^. Simultaneous seven-color GROND 
photometry further probes the heating of this companion, and shows the presence of a flaring infrared excess. 
We have modeled the quiescent light curve, constraining the orbital inclination and masses. Simple heated 
light-curve fits give il^NS = 2.7 M0, but show systematic light-curve differences. Adding extra components 
allows a larger mass range to be fit, but all viable solutions have Mns > 2.1 Mq. If confirmed, such a large 
Mns substantially constrains the equation of state of matter at supernuclear densities. 
Subject headings: gamma rays: stars — pulsars: general 



1. INTRODUCTION 

The bright 7-ray source 2FGL J131 1.7—3429 has been 
known, and unidentified, sin ce the early E GRET mission 
(Fichtel et al. 1994). After Romanil (120121 hereafter R12) 
discovered an optical counterpart with a = 5626.0 s 
(93.8 min) orbital period and evidence for very strong pul- 
sar heating, iPletsch et al.l (12012) managed to discover 2.5 ms 
Doppler-mo dulated pulsations in the Fermi Large Area Tele- 
scope (LAT. lAtwood et aLl 1200 9*) 7-ray photons. Thus with 
orbital constraints, 7-only millisecond pulsars can be detected 
by the LAT, although PSR J1311-3430 was later found to 
also be intermittently visible in the radio (Ray et al. 2012, 
in prep.). Additional optical and X-ray observations sup- 
porting the orbital v ariability found above are described by 
lKataokaetai](l20T2h . 

The 7-ray pulsar ephemeris gives an orbit of ans sin i — 
0.0105811t-s for a minimum companion mass of ^ 8 x 
1O~^M0. This measurement alone does not significantly 
constrain the neutron star mass, which is of particular inter- 
est, since Ivan Kerkwiik. Breton. & Kulkarnil (1201 IL hereafter 
vKBK) found evidence that the original, Pt, = 9.2 hr "black 
widow" (BW) binary PSR B1957H-20 may host a remarkably 
heavy neutron star- with A/ns = 2.4±O.12M0. If BW pulsars 
as a class are massive, PSR Jl 131—3430, with its ultra-short 
orbital period (the smallest of any rotation-powered pulsar), 
provides interesting opportunities for a mass determination. 

We have been able to obtain exploratory spectroscopy and 
additional photometry that constrain properties of the PSR 
J131 1—3430 system. These data show that the companion 
is a bloated, Roche-lobe filling substellar object whose pho- 
tosphere is He dominated, with no detected H. While the neu- 
tron star mass estimate is quite model dependent, the data al- 
ready imply that it is quite large. We discuss here the obser- 
vations, the implications, and the prospects for further refine- 
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ment of the system parameters. 

2. KECK LRIS SPECTROSCOPY 

Six consecutive 300s exposures of J1311— 3430 were ob- 
tained with the Low Resolution Imaging Spectrometer (LRIS; 
Oke et al. 1995) at the Keck-1 10 m telescope on 2012 May 17 
(UT dates are used herein; MJD 56064). We used the 1" long 
slit at the parallactic angle (Filippenko 1982) and the 5600 A 
dichroic splitter. In the blue camera the 600/4000 grism pro- 
vided coverage to < 3400 A at ~ 4 A resolution; the red cam- 
era employed the 400/8500 grating for ~ 7 A resolution. The 
seeing was good 0.8") and conditions were nearly pho- 
tometric. Standard processing and optimal extraction were 
applied; fluxes were calibrated against the red and blue stan- 
dards BDh-26°2606 and BDh-28°421 1, respectively. 

These observations covered the phase of maximum light 
(0.7 < cj)B < 1.1), so the spectrum is very blue (Fig- 
ure 1). Although the temperature inferred from the color 
is >10,000K (see §3), no Balmer lines are visible, im- 
plying an unusual surface composition. This is in con- 
trast to the BW companio n s of PSR B1957H-20 (vKBK; 
lAldcroft, Romani, & Cordii [T992b and 2FG J2339-0533 
(Romani & Shaw 20 fl]), where relatively normal companion 
spectra are seen, with spectral classes appropriate to their 
(cooler) surface temperatures. After shifting the individual 
spectra to the rest frame using the sinusoidal approximate 
Doppler curve we obtain a phase-averaged spectrum domi- 
nated by narrow neutral He lines. We infer that this is the 
H-stripped remnant of an evolved stellar core. 

The lower panel of Figure 1 shows the blue region of 
the normalized spectrum compa red with a He-dominated 
model from the atmosphere grid of lJeffery. Woolf. & Pollacc"ol 
(12001b . smoothed to the resolution of the Keck data. The cor- 
respondence is excellent, with the exception that the Balmer 
lines, well visible in the model despite the 10~^ H number 
abundance, are undetectable in PSR J131 1—3430. This im- 
plies a maximum H abundance n(H) < 10~^. 

We have used the model grid to estimate the effective tem- 
perature and gravity of this phase-averaged spectrum. The fit 
minima give Tcff = 12, 000 ± 1000 K and logg = 4.6 ± 0.2. 
Measurements of the equivalent width (EW) of individual 
lines are consistent with these conclusions; for example, the 
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Fig. 1 . — Combined, phased Keck LRIS spectram. Upper Panel: Full spec- 
trum. Lower panel: Spectrum showing He I line domination and the good 
match to the JefferY^Woolf. & Pollacco (2001) model atmosphere (blue). 

ratio EW447i/liiW448i = i.^ ± u:z and the absence of He II 
support the T^ff estimate, while the ratio EW4471/EW3819 = 
1.97 ± 0.2 accords with the log g value. A few species, such 
as Si A4552 and Ni A3995, suggest some flux at lower val- 
ues of log g. This is not very surprising, since the compan- 
ion is likely near-Roche-lobe filling with appreciable surface- 
gravity variation. Higher signal-to-noise ratio (S/N) spectra 
and synthetic composite spectral models will be needed to 
study the atmosphere's surface variation. 

The highly unusual He-dominated, low \ogg spectrum 
compromised attempts to find suitable cross-correlation tem- 
plates. Accordingly, we adopt the best-fit atmosphere model 
to measure the radial velocities, using the RVSAO package 
in IRAF. All observed comparison stars were poor spectral 
matches and gave significantly lower correlation coefficients. 
However, all templates gave very similar radial velocities. 

With the phase fixed by the pulsar-determined orbital 
ephemeris, we can fit a simple sinusoid to these velocities, 
obtaining amplitude Kohs = 609.5 ±7.5 km s^^ and average 
radial velocity T = 62.5 ± 4.5kms"\ with Xu 1-3- At 
face value this gives a mass function 

/2,obs = = ^ = 1.54 ± 0.06 Mo 

and a very large mass ratio gobs = (-K^obsZ-fi'Ns)'^ — 172.1 ± 
2.1. However, the observed spectrum tracks the center of 
light (CoL) from the heated face of the companion, which has 
a smaller radial-velocity amplitude than the companion cen- 
ter of mass (CoM). For PSR B1957H-20, vKBK estimate the 
correction factor as A'cor = 1-09 (i.e., i^coM = A'cor^obs)- 
If we adopt this value of K^oy and the inclina tion constraint 
i < 8 5° from the lack of X-ray eclipses (R12: iKataoka et aP 
I2OI2I) . we derive A/ns > 2.03 Mq. This is an interestingly 
high mass for equation-of-state (EoS) constraints. 

Ho wever, Kco i depends on the heating across the surface. 
Also, iRomani & Shaw. (.2011.) find that for the similar BW- 
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type binary J2339— 0533, the heating introduces a significant 
non-sinusoidal component to the CoL curve, and it should be 
fit. Any further restrictions on the inclination i greatly narrow 
the pulsar mass range. These factors can be probed with light- 
curve modeling, so we have sought to improve the photometry 
to allow better constraints on the system masses. 

3. GROND PHOTOMETRY 

R12 measured the companion light curve from WIYN, 
SOAR, and archival VLT 3-band (g'r'i' or BVI) photome- 
try. However, the colors were not simultaneous; the largest 
dataset (SOAR Optical Imager g'r'i') covered each filter on a 
separate night. These data showed significant epoch-to-epoch 
variability and large (up to 4 mag) flares, so precise measure- 
ments of the instantaneous colors were not available. Ac- 
cordingly, we targeted J 1 3 1 1 -3430 with the GROND system 
on th e MPI/ESO 2.2m telescope on La Silla (GreinereLal] 
|2008| ), obtaining simultaneous g'r'i'z'JHK images starting 
2012 July 09. The observations lasted 2.5 hr 1.5 orbits) 
from MJD 56117.956 to 56118.059, covering two maxima 
(pulsar inferior conjunction). Twelve observation blocks were 
made, each with 4 x 115 s dithered optical pointings and 
8 X 60 s near-IR exposures. The seeing was ^ 1" and the air- 
mass increased slowly. Near maximum light, individual g'r'i' 
exposures had adequate S/N, allowing a temporal resolution 
of ^ 195 s, while at minimum and for all z' images 4 dithered 
exposures were combined to achieve adequate S/N. To ensure 
the best possible calibration of the optical bands, the target 
was revisited on 2012 August 10 with SDSS field exposures 
bracketing the J131 1—3430 observations, allowing an accu- 
rate calibration of the field and target magnitudes; the calibra- 
tion systematics are small compared to the photometric errors, 
and are included in the error flags of Figure 2 {ag' « 0.02, 
(Tr' « 0.03, ai' « 0.05 mag at maximum). 

J131 1—3430 had relatively low near-IR S/N and was de- 
tected at high significance only near maximum light; we thus 
made photometric measurements on an image stacked from 
frames covering phase ±0.17. The near-IR data were cal- 
ibrated using 2MASS stars in the field, and transformed to 
AB magnitudes for spectral energy distribution (SED) fitting 
together with the optical data. These peak magnitudes were 
J = 20.6 ± 0.1, H = 21.1 ± 0.3, and K = 20.9 ± 0.3. All 
magnitudes were corrected for the filter-specific dust-map- 
estimated Galactic extinction {Ay = 0.173 mag), using the 
Schlafly & Finkbeiner (201 1) calibration. 

Figure 2 summarizes the GROND photometry. The left 
panel shows the g'r'i' z' light curves. The lower portion dis- 
plays the g' — i' color, blue at maximum and redder at min- 
imum. The upper right shows the optical through near-IR 
SED at maximum light (—0.1 < (j>B < 0.1) and the optical 
SED near quadrature (1^)3 = 0.2-0.3, 0.7-0.8). We see that 
at maximum the colors are comparable to those of a B7-B8 
star in the blue, but there is a large near-IR excess. Although 
the pulsar heating makes the star multi-temperature, this very 
large excess suggests a large emitting area at low tempera- 
ture. One source larger than the companion photosphere is 
the evaporative wind, where pulsar power may be reprocessed 
into the optical-IR. The optical light curves show appreciable 
fluctuations on short timescales. The lower-right panel sum- 
marizes this variability with the g' — i' color plotted against 
the g' magnitude. Again the star is redder at fainter magni- 
tudes. Interestingly, if we flag points that are brighter than the 
"quiescent" light-curve flux (see below), most of these points 
shift (green lines) to larger g' — i' (redder) relative to the ex- 
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Fig. 2. — GROND photometry. Left: g'r'i' z' light curves, showing red- 
dening during minimum (pulsar superior conjunction). Upper right: The 
companion SED at maximum (phases —0.1 to 0.1. blue points) and quadra- 
ture (phases 0.2-0.3 and 0.7-0.8, red points). Comparison with a B8 star 
shows a large IR excess. Bottom right: J1311— 3430 color variations. The 
redder spectrum at minimum is visible. Green lines and red triangles mark 
when an epoch has Ag' > 0.2 mag over the quiescent magnitude. Such 
flares show substantial increases in g' — i'. 

pected quiescent color. An extreme example is the very large 
i' flare recorded by R12. We conclude that the principal light- 
curve fluctuations have a red SED. This implicates modula- 
tion in reprocessed pulsar emission by the variable wind off 
the companion. An alternative fluctuation site, stellar flares 
on the tidally locked companion, would in contrast be very 
hot (blue). 

4. ELC MODELING AND PARAMETER ESTIMATION 

The GROND data provide an instantaneous color reference 
for the longer, higher precision, but non-simultaneous SOAR 
Hght curves from MJD 56008.2-56010.4. Just as the GROND 
data exhibit red flares above the "quiescent" flux, comparing 
SOAR g'r' photometry from adjacent orbits reveals several 
~ 20min flaring periods. In Figure 3 we plot the phased 
SOAR/GROND photometric points. During the first period 
all g'r' measurements and errors are shown — "flaring" points 
(at0B — 0.15-0.35 and 0.45-0.8) are plotted as crosses. Dur- 
ing the second period we display only the quiescent points, 
but include the GROND i'z' quiescent data. All SOAR i' 
were affected by strong flaring. 

To constrain the system properties we have fit these "qui- 
escent" light curves along with the radial-velocity measure- 
ments, using the "E clipsing Light Curve" (ELC) code of 
lOrosz & Hauschildl 0200 0). Typically, the code computes 
filter colors from atm ospheres of the "NextGen" library 
(iHauschildt et al.ll999h . but we found that these did not match 
the data well, plausibly due to the He domination of the pho- 
tosphere. Instead, we fit using blackbody emissivities which 
gave reasonable light curves and colors. All fits used the 7- 
ray-determined phasing and pulsar projected semimajor axis. 

There are five basic fit parameters: the system mass ra- 
tio q = M^a,/Mc, the pulsar irradiation modeled here as 



an isotropic flux Lx, the companion's underlying tempera- 
ture Ti and Roche-lobe fill factor /i, and the orbital incli- 
nation i. Lx is quite well constrained by the g' — r' and 
r' — i' colors near maximum light, giving log (Lx) = 35.3. 
This is in accord with the spectroscopically estimated Toff = 
1 2,000 K near maximum. We also find that the secondary is 
very close to Roche-lobe filling in all acceptable fits, so we 
set /i = 0.99. We therefore adjusted the three remaining 
parameters to fit the light curves and spectral points. The 
very small photometric errors near maximum, coupled with 
stochastic variability, meant that even the best fits had formal 
X^/dof f« 9, so the Ti and i error estimates come from the 
range giving a Ax^ increase of (x^/do/)min around the fit 
minimum. For q (and A/ns) the value (for a given fit i) is 
controlled by the radial-velocity measurements. These had 
(x^/(io/)inin w 1.5, SO the error range was estimated from 
this smaller xl increase. 

For both PSR B 1957-H20 (iRevnolds et al.1 l2007h and 
J2339-0533 ( lRomani & Shawl 1201 Ih the photometric data 
are limited, and simple models with a pulsar-heated hemi- 
sphere gave adequate light-curve fits. In contrast, the flat 
light-curve maxima of J131 1—3430 produce relatively poor 
fits to such models (R12). The best fit (Table 1) has a 60° in- 
clination, resulting in a rather large pulsar mass, but the peak 
is too narrow, the minimum is too bright by ^ 1 mag, and the 
predicted /fobs is 2.2a below the measured value. 

Adding equatorial hot spots to the companion surface 
bracketing the LI point (30° radius at 35° and —65°) can 
broaden the light-curve peak. This mimics equatorially con- 
centrated heating, which may be plausibly invoked, as there is 
strong equatorial concentration in pulsar wind nebula (PWN) 
flows, giving rise to PWN tori. Alternatively, these asym- 
metric components may represent reprocessed pulsar flux in 
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Fig. 3. — Top: SOAR (filled points) and GROND (open points) photom- 
etry; g'r' "flare" measurements are flagged as crosses. Points during the 
second period exclude flare epochs. Curves are from the "LI Spot" ELC 
fit. Bottom: Keck radial-velocity measurements, aligned with the light curve. 
Residuals are shown amplified by a factor of 10. The curve is from the same 
"LI Spof fit. 
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TABLE 1 
ELC Fit Parameters 



Parameter 


Basic LC 


LI Cold Spot 


Eq Hot Spots 




60.4±0.4 


67.3±0.3 


57.9±0.3 


Ti[K] 


3440±50 


<2000 


< 1600 


<? 


179.7±3.9 


177.1±3.2 


180.2±3.3 


Mns[M0] 


2.68±0.14 


2.15±0.11 


2.92±0.16 


^cor 


1.06 


1.04 


1.06 



a companion wind outflow. With these added components 
the model is a good match to the g'r'i' light curves. The 
heated equator allows slightly lower log(Lx) — 35.0, but 
the KcoY = 1-06 is large and the best-fit inclination is small, 
giving a very (possibly unphysically) large mass. 

If, in contrast, we cool the region near LI, by applying a 
40° spot with 1/2 the local T, we also flatten the maxima. In 
this case, the CoL moves toward the CoM, resulting in a small 
KcoT = 1-04 and larger i. These two effects allow a neutron 
star mass of 2.15 Mq. The fit quality is good, albeit with 
higher than the equatorial spot model. This scenario gives 
the smallest A/ns of any viable fit. Such decreased LI flux 
might plausibly arise from large limb and gravity darkening 
effects in the He-dominated atmosphere. We thus attempted 
to flatten the light curve by invoking extreme gravity darken- 
ing coeffic ients, as have been claimed f or some semidetached 
binaries dNakamura & Kitamural[T992l) . This allowed inter- 
mediate inclination and masses; for example, for Tgi = 0.5, 
we find i = 63° and A/ns = 2.30 M0. For all cases the fit 
models underpredict the z' flux. This is Hkely related to the 
extra IR component seen in Figure 2. 

5. DISCUSSION AND CONCLUSIONS 

It is clear that additional physics is needed to fully model 
the light curve. Certainly, specific intensities from an ap- 
propriate He model grid would be useful. However, a good 
fit may also require some asymmetry in the heating, or re- 
processing of pulsar flux to produce optical/IR emission from 
the companion wind. We conclude that, at present, model 
assumption systematics dominate the statistical fit errors and 
preclude accurate mass determination. Nevertheless, the 
large mass function of PSR J131 1—3430 virtually guaran- 
tees a large Mns- In this it joins B1957H-20, supporting the 
suggestion of vKBK that BW pulsars as a class may have 
high masses. The challenge is that secondary spectroscopy- 
derived estimates for these BW masses inevitably have sub- 
stantial systematic uncertainty. In this sense, the mass esti- 
mates cannot match the "gold standard " timing-based Mnr = 
1.97± 0.04 Mq for PSR J1614-2230 (iDemorest et al.l2010h . 
However, since EoS constraints tighten rapidly as the mini- 
mum required A/ns exceeds 2 Mq, it is worth considering the 
prospects for refining the J131 1—3430 mass estimate. 

Figure 4 shows the present situation in the mass-mass plane. 
Compared to PSR B1957H-20, J1311-3430 has a larger ob- 
served mass function (1.54 M0 vs. 1.34 M0). However, our 
models give substantially smaller Kcoi than assumed for PSR 
B1957H-20, so in this sense our mass fits are conservative. The 
remaining systematic uncertainties in the light-curve shape al- 
low inclinations 57° < i < 67°. This large range dominates 
the systematic mass uncertainty. The best hope for eliminat- 
ing otherwise plausible models and reducing the systematics 
lies with additional spectroscopy. Full-orbit coverage will, 
at minimum, substantially decrease the uncertainty in q. We 
also expect to detect the non-sinusoidal radial-velocity com- 
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Fig. 4. — Mass constraints on the PSR J131 1—3430 system. The pulsar 
mass function gives the solid diagonal lines for vaiious inclinations i. The 
dotted vertical line gives the minimum from the companion mass function. 
The Kcovv increase of the observed radial velocity is critical: we show the 
physically allowed range from the Keck data (dotted diagonal red lines, la- 
beled by the appropriate value of q for Jl 31 1—3430), along with the range 
for PSR B 1957+20 preferred by vKBK(_ft'corr = 1.06-1.13; shaded region). 
Fit regions for the three models in Table 1 are indicated with blue ellipses. 

ponent, and with sufficient spectral resolution, the variations 
in absorption-line profile expected as the visible heated sur- 
face varies. Additional simultaneous multicolor photometry 
from many epochs could further isolate the photospheric light 
curve from a variable wind component. Deep multicolor im- 
ages can constrain the residual illumination at pulsar supe- 
rior conjunction, with additional i constraints. Of course, im- 
proved models are needed to fully exploit such data. 

Our data also illuminate several other aspects of this unique 
system. The IR-dominated variability likely probes the ex- 
tremely dense wind flux. Correlation with periods of radio- 
pulse visibility could be very revealing. The ultra-low-mass 
He companion points to unusual evolution, with substan- 
tial mass transfer and irradiative stripping of the core of an 
evolved secondary. Benvenuto, De Vito, & Horvath (2012|) 
describe such a scenario, giving parameters remarkably simi- 
lar to those of J131 1—3430. This model's slow mass transfer 
should allow substantial neutron star mass growth. Interest- 
ingly, for the large neutron star masses indicated here, the EoS 
must be very "stiff." This means that for A/ns ~ 2.5 Mq, we 
can hav e a moment of inertia as large as ~ 4 x lO''^ g cm^ 
jLattimer & Prakash ,200 7). This increases the spindown 
power of J1311-3430 from £; = mtl = 5 x lO^^ergs"! 
(for I45 = 1), accommodating the large heating flux inferred 
from the peak colors and spectroscopic Toff, even without 
equatorial concentration. 

In summary, PSR J131 1—3430 presents a number of 
unique properties that make it a potential Rosetta Stone for 
irradiation-driven evolution in tight binaries. If we can reduce 
the allowed model space, it seems likely to provide some of 
the most important constraints on the EoS of dense matter 
And, even if PSR J 1 3 1 1 -3430 or PSR B 1957H-20 do not pro- 
vide the final word on the EoS, it is worth remembering that 
Fermi has proven to be an excellent finder of short- Pb BW 
systems, including detection via 7-ray-blind searches. Future 
studies of this extreme population may well provide a bullet- 
proof case for a stiff EoS. 
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